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ABSTRAGT 

Lyct photons scattered by neutral hydrogen atoms in the circumgalactic media or produced in the 
halos of star-forming galaxies are expected to lead to extended Lyct emission around galaxies. Such 
low surface brightness Lyct halos (LAHs) have been detected by stacking Lyct images of high-redshift 
star-forming galaxies. We study the origin of LAHs by performing radiative transfer modeling of 
nine ^ = 3.1 Lyman-Alpha Emitters (LAEs) in a high resolution hydrodynamic cosmological galaxy 
formation simulation. We develop a method of computing the mean Lya surface brightness profile of 
each LAE by effectively integrating over many different observing directions. Without adjusting any 
parameters, our model yields an average Lya surface brightness profile in remarkable agreement with 
observations. We find that observed LAHs cannot be accounted for solely by photons originating from 
the central LAE and scattered to large radii by hydrogen atoms in the circumgalactic gas. Instead, Lya 
emission from regions in the outer halo is primarily responsible for producing the extended LAHs seen 
in observations, which potentially includes both star-forming and cooling radiation. With the limit on 
the star formation contribution set by the ultra-violet (UV) halo measurement, we find that cooling 
radiation can play an important role in forming the extended LAHs. We discuss the implications and 
caveats of such a picture. 

Subject headings: cosmology: observations — galaxies: high-redshift — radiative transfer — scattering 
— intergalactic medium 


I. INTRODUCTION 


The Lya line is an important cosmological tool for 
studying star-forming galaxies in the young universe, and 
has been found to aid in the detection of high redshift 

Ouchi et al.||2QQ8 Gawiser 


galaxies dRhoads et al.||2QQ3l \0 
et al.||2QQT " Guaita et ar][2in(^ 


As ionizing photons are 
emitted from young stars, they ionize neutral hydrogen 
in the surrounding interstellar medium, and are likely to 
be re-emitted as Lya phot ons following recombination 


(Partridge & Peebles 1967). After they escape the in- 
terstellar medium surrounding their parent stars, they 
are predicted to undergo resonant scattering with neu¬ 
tral hydrogen gas in the surrounding medium as a result 
of the radiative transfer process, diffusing ou t both spa¬ 


tially and in frequency (Zheng et al. 2011a). As such, 
extended halos of neutral hydrogen around these Lya 
emitters (LAEs) are predicted to be illuminated by scat¬ 
tered Lya photons. 

Many theoretical studies have predicted the existence 


shift galaxies (e.g., Tasitsiomi||2006| 

Laursen & Sommer- 

Larsen||2007 Laursen et al.||20U9| L 

bjkstra V Loeb||2UU9t 

Barnes & Haehnelt|20101|Barnes et 

al.|2011t|Zheng et al.l 

2010[ 201 la|). While these LAHs are predicted to gener- 


z >2^ their presence can be revealed by stacking tens to 
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hundreds of narrow band images of high redshift LAEs 


(iFynbo et al.|2Q Q3 S teidel et ai.|2QlT Zheng et al. 2011a 

Matsuda et al.pff L^ llVIomose et al.|2(J14p . Although ob- 
here s 


servationally there seems to be a consensus in favor of 
their exis tence, there are also re ports of null detections 
of LAHs. Feldmeier et al. (2013) find marginal evidence 
and no evidence of LAHs for 2 : ^ 3.1 and z ^ 2.1 LAEs, 


while Jiang et al. (2013) find evidence of LAHs based 
on results using stacked images of LAEs at redshifts of 
5.7 and 6.6. Such contradicting results may be ca used 
by small number statistics or unknown systematics (Mo- 

et al.llMl . 


mose < 


T'he shape and size of LAHs can yield insights into the 
spatial distribution and kinematic properties of the cir- 
cumgalactic and int ergalactic medium surrounding LAEs 
(jZheng et al.|2011a|). The shape can also be used to con¬ 
strain cosmic reionization, with rei onization leading to 
steep er surface brightness profiles (Jeeson-Daniel et al. 
2012). Detailed theoretical studies of LAHs can help in 
understanding their origin and properties. 

The aim of this pa per is to apply a Monte Carlo ra- 
diative transfer code (Zheng & Miralda-Escude||2002|) to 
study diffuse Lya halos surrounding 2 : ^ 3.1 star-forniing 
galaxies in a high-resolution galaxy formation simulation. 
By comparing with observational data, we hope to gain 
insight about the origin and composition of these diffuse 
LAHs. 

This paper is divided into several sections. In Sec¬ 
tion 2, we describe the modeling method and the con¬ 
struction of the average Lya surface brightness profile 
for each model LAE. Our main analyses and results are 
presented in Section 3, with comparisons to observations 
and a discussion of possible constraints imposed by the 
profile in the UV band. Finally, we summarize our re¬ 
sults and discuss the implications in Section 4. 
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Fig. 1.— Neutral hydrogen column density maps for the model LAEs in our analysis. Each image is 224 kpc (physical) on a side. The 
column density is computed by integrating over the whole box along the line of sight (224 kpc physical). The black contours are drawn 
at 10^®-^cm“^, within which are regions corresponding to DLAs. The white contours are drawn at lO^^-^cm”^. Extended filamentary 
structures of neutral hydrogen are seen, which are connected to the extended Lya emission discussed in this paper. 


2. MODEL 

2.1. Lya Radiative Transfer Modeling of Simulated 
Star-forming Galaxies 

Our Lya radiative transfer modeling of simulated star¬ 
forming galaxies is based on a cosmological simulation 
with the adaptiv e mesh refinement (AMR ) Eulerian hy- 
dro code Enzo dBrya n fc NormanI |2QQQ| Ijoung et al 


2009 


as detailed^ in |Cen 
In brief, a region o; 

^3 • 


(|2U12|| 

4 conic 


and [^n ^ 

comoving size 21 


Zheng 


20h“^Mpc'^ in a low-resolution simulation (with a box 
size of 120h“^Mpc comoving on a side) is chosen to be 
resimulated at high resolution. The resimulation has a 
dark matter particle mass of 1.3 x 10^h“^M q and the 
mesh refinement ensures a spatial resolution better than 
lllh“^pc (physical). The resimulation includes an ion¬ 
izing UV background and the self-shielding of the gas, 
metallicity-dependent radiative cooling, molecular hy¬ 
drogen formation, star formation, and supernova feed¬ 
back. The mass of a star particle is typically ^ IO^Mq. 
The simulation assumes a spatially flat ACDM model 
with the following cosmological parameters: Qrn = 0.28, 
fib = 0.046, Ho = 100hkms“^Mpc“^ with h = 0.70, 
cTg = 0.82, and Ug = 0.96. 

The simulation has been used to study the kinematic 
properties traced by unsaturated metal lines in damped 
Lya s ystems, wh ich is in good agreement with observa¬ 
tions ( Cen||2012 ). The simulation has also been applied 


to study the partition of stell ar light into optical a nd 


infrared light (|Cen||2011|). In |Cen & Zheng (|2013j), 
model of Lya blobs (LABs) is developed based on L 


model of Lya blobs (LABs) is developed based on Lya 
radiative transfer modeling of the simulated star-forming 
galaxies in massive halos, and the observed relation be¬ 
tween Lya luminosity and LAB size and LAB luminosity 
function at 2 : ^ 3.1 have been successfully reproduced. 
In this paper, we select from the simulation 9 2 : = 3.1 
star-forming galaxies in halos of mass to study 

the properties of LAHs associated with them. The value 
of is chosen as a starting point for our analysis, 

and is within current constraints on LAE halo mass of 


10^^^^Mq ( Ouchi et al.|2QlQ ). The mean stellar mass of 


these nine galaxies is about 2.9 x IO^^Mq. _ 

We implement a Monte Carlo code developed by [Zheng] 


& Miralda-Escude (2002) for the Lya radiative transfer 
calculation in extended neutral hydrogen distributions 


surrounding our simulated LAEs. Thi s code has been ap' 
plied to study LAEs and LABs (e.g., Zheng et al. 1 12010^ 
2011b|a Cen & Zheng 2013 Zheng If: Wallacej |2014p ' 


For each galaxy, we store the relevant quantities from 
the simulation in a uniform cubic grid of 4i?vir on a side, 
with cell size 319 pc (physical, corresponding to 0.04"). 
Here i?vir is the virial radius of the host halo, which is on 
average ^ 56 kpc (physical) for the nine halos 

we consider. The quantities include the Lya luminosity, 
neutral hydrogen density, temperature, and velocity. The 
Lya luminosity is separated into star-formation and cool- 
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Fig. 2.— Lya surface brightness images for all nine model LAEs in our analysis. Each image is 224 kpc (physical) on a side and has 
been smoothed by a 2D Gaussian kernel with a FWHM of 1.32'' to match the observation setup. Isophotal contours are drawn in black 
at the limit of observational detection, 10“^^ergs“^ c m~^ arcsec~^. The da shed circle in each panel has a radius of 40 kpc, beyond which 
the systematic effects in the image stacking analysis in|Momose et al.|(|2014|) start to become important. The images show a rich diversity 
in structure at surface brightness levels below the detection tJiresJiolcl, wJiicn contributes to the extended LAHs. 

ing contributions. The Lya luminosity from star forma- 
tion is computed as = lO^^[SFR/(M 0 yr“^)]erg s 


(Furlanetto et al. 2005), where SFR is the star forma- 


as viewed alo ng the chosen direction (Zheng & Miralda- 


Escude 

2002 ), 

Final. 

ly, we 


tion rate in the cell. The Lya luminosity from cooling 
radiation is computed from the de-excitation rate, which 
depends on neutral hydrogen density and temperature 
that are computed self-consistently by following the rel¬ 
evant species in a non-equilibrium fashion. 

Each photon launched from a cell is assigned a weight, 
calculated by dividing the total Lya luminosity of the 
cell by the number of simulation photons lauched from it. 
Such cell-dependent weights are accounted for in comput¬ 
ing the hya surface brightness profiles. The scatterings 
of the photon with neutral hydrogen atoms on its way out 
and the corresponding changes in position, direction, and 
frequency are tracked until it escapes the grid boundary. 
We record the initial position of each photon, the posi¬ 
tion of the last scattering, the direction and frequency 
after the last scattering, and the fractional contribution 
of cooling radiation to its total luminosity. This informa¬ 
tion is used to compute a mean surface brightness profile 
for each LAE, averaged over all directions (see § 2 . 2 ). 

At each scattering, we also compute the escape prob¬ 
ability towards a fixed direction and collect the es¬ 
caping hya photons onto an integral-field-unit-like 3- 
dimensional array with pixel size the same as the cell size, 
which allows us to construct a Lya image of each LAE 


tic medium (IGM) outside of the box and the interstel¬ 
lar medium (iSM) in star-forming regions on Lya emis¬ 
sion, f ollowing the approximate methods in |Cen fc Zheng] 


(2013). 

in brief, for each photon escaping the box at frequency 
z/, we calculate the scattering optical depth from the 
edge of the box to an observer at 2 ; = 0 using the redshift- 
dependent IGM hydrogen density and apply a factor of 
correction for the IGM absorption. While such 
a correction neglects the differences in the IGM along 
different directions, it serves our purpose of introduc¬ 
ing an overall average effect of the IGM. We also ap¬ 
ply an effective ISM dust attenuation to the intrinsic 
Lya emission by multiplying the luminosity represented 
by each simulation photon by a simple e~^ factor, with 
r = O. 2 [SER/(M 0 yr“^)]^-^. This is loosely motivated 
by the observational tren ds of higher dust attenuation in 
galaxies wi th higher SER ( [Brinchmann et al.|2QQ4 Zahid 


et al.|2Q13D , and the power-law index follows the slope of 


the metal column density dependence on SER in the sim¬ 
ulation. The dust extinction can be thought as applied 
to the ionizing photons around the H II regions, lowering 
the luminosity of Lya emission coverted from ionizing 
photons through recombination. The factor is also in- 
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tended to absorb uncertainties in the galaxy formation 
simulation (e.g., in the predicted SFR). Our methods 
of applying IGM and ISM abso rption are the same as 
adopted in Cen & Zheng (2013), in which the LAB lu¬ 
minosity function and lurninosity-size relation have been 
successfully reproduced, suggesting that the approximate 
treatments work well in capturing the major IGM and 
ISM effects and in absorbing model uncertainties. In our 
current work of LAHs, we do not adjust any parame¬ 
ters and simply use the direct outputs of the radiative 
transfer model to compare to observations. 

We show in Figure the column density distribution 
of neutral hydrogen gas around each of our model LAEs, 
viewed from a fixed direction. The black contour curves 
are drawn at lO^^’^cm”^, delineating regions correspond¬ 
ing to damped Lya systems (DLAs). DLAs represent re¬ 
gions extremely opaque to Lya photons, while Lya pho¬ 
tons can be significantly scattered in regions of much 
lower column densities (e.g., above 10^^cm“^). The im¬ 
ages reveal extended, filamentary structures of neutral 
hydrogen connecting regions of high column densities. 
Scatterings of Lya photons off hydrogen atoms in these 
structures leave signatures in the resulting Lya surface 
brightness distributions. 

Figure shows the corresponding Lya images of the 
nine LAEs in our analysis. The isophotal contours in 
each image correspond to I0“^^ergs“^ cm“^ arcsec“^, 
about the surface bright ness threshold for detecting in¬ 
dividual 2 ; ^ 3.1 LAEs (jOuchi et al.||2008|. These im¬ 
ages reveal a rich degree of structures and a variety of 
morphologies at fainter surface brightness levels, which 
allows the LAHs to be revealed by the stacking analy¬ 
sis. The surface b rightness distribu tion depends on the 
viewi ng angle (e.g. [Zheng e t al.|2QlQ|). The stacked image 
(as in IVlomose et al.||2(J14|) for a sample of LAEs comes 
from averaging many images from galaxies of random 
orientations. Eor the relatively small number of galaxies 
modeled here, we can form the mean surface brightness 
distribution by viewing each galaxy from many different 
observing directions. 


2.2. Computing the Mean Lya Surfaee Brightness 
Profiles of Model LAEs 

Eirst, let us consider an observer located along a di¬ 
rection k from one of our model galaxies. The aver¬ 
age surface brightness SB(i?, k) at a projected radius R 
(physical) to the galaxy center as seen by the observer 
can be computed as (assuming a spatially flat universe) 


SB(R, k) 


AL{lPz)-^ 

mwHaoTWMV 


( 1 ) 


where ALl is a small solid angle centered around k, A = 
2t:RAR is the area of a small annulus around R, and AL 
is the total luminosity of escaped Lya photons falling into 
Aft and with projected last scattering position within A. 
The quantity AL{1 z)~^/{AQD^) is the corresponding 
flux (erg s“^ cm“^) the observer at a comoving distance 
Dc receives, with the two factors of I + 2 ; from energy 
redshift and time dilation. The quantity A(1 zfi^ jD^ 
is the solid angle extended by the annulus, seen by the 
observer, with the (I+ 2 :)^ factor converting physical area 
to comoving. 

The surface brightness profile at projected radius R 


averaged over all observing directions is then 


(SB(i?)) 


1 

47r 


/ 


SB(R, k)dn. 


( 2 ) 


Given Equation ([B, for the annulus at a given R and AR, 
we only need to obtain the average of AL/ALL over all 
observing directions for computing the integral in Equa¬ 
tion (©■ Denoting the total Lya luminosity from this 
annulus as La^ we then have La = / AL/ALtdLt. In the 
limit of an infinite number of observing directions, we 
obtain 


(SB(i?)) = 


La 

47rA(I -1- zY 


(3) 


This equation is the basis of computing the mean Lya 
surface brightness profile for each model LAE, along 
with the information we record for the escaping Lya 
photons. Eor the mean Lya surface brightness at pro¬ 
jected radius R, instead of producing images viewed from 
many directions, we only need to obtain the sum of the 
total Lya luminosity La for photons whose projected 
radii fall into the annulus around R {R A AR/2 with 
an area A). The projected radius R^ of a photon is 
computed from its escaping direction k^ and its position 
of last scattering rig (with respect to the galaxy center) 
as Rry = — (k^ • rig)^. We have verified that the 

method gives the same results as that from averaging 
images over many observational directions. 

Our simulation has a much high er resolution (cell siz e 
of ^0.04") than the observation in Momose et ah (|2QI4 ). 
To m i mic th e smoothing effect in the images of|lVlomose 


et al. (2014), we obtain the final surface brightness prohle 


by convolving the resulting Lya image with a 2D Gaus¬ 
sian kernel with a full width at half maximum (EWHM) 
of 1.32", corresponding to 10.3 kpc(physical) at 2 : = 3.1. 


3. RESULTS 

We first present the results on the mean Lya surface 
brightness profile in our model. We then decompose the 
mean profile in various ways to study its origin. Einally, 
we compare our UV profile to observations in order to 
further constrain the relative contributions of cooling and 
star-forming emission. 


3.1. The Mean Lya Surfaee Brightness Profile 

The left panel of Eigure shows the mean surface 
brightness profiles from the smoothed images of the nine 
model LAEs in our analysis. The overall mean of the 
nine profiles is plotted in black. Eor each LAE, the 
central profile (e.g., R < 15 kpc) is largely determined 
by the smoothing kernel (point spread function; PSE). 
The peak surface brightness at the center shows a small 
variation among the nine individual profiles, around (1- 
3)xl0“^^ergs“^ cm“^ arcsec“^. 

At large projected radii, each LAE shows an extended 
profile, which can be identified as the diffuse LAH. The 
profile at large radii is much flatter than the central part. 
The surface brightness level of this extended component 
displays a substantial variation among the nine individ¬ 
ual LAEs, as large as two orders of magnitude. 

In the right panel of Eigure we compare the overall 
mean pr ofile of the nine mod A LAEs with the one de¬ 
rived by Momose et al. ( 2Q14| ) from stacking the Lya 
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Fig. 3. — Lya surface brightness profiles for the model LAEs in our analysis. Left: Lya surface brightness profiles for individual model 
LAEs, with the average over all nine LAEs shown as the black curve. While the nine LAEs have similar surface brightness levels at small 
radii, they display large variations at large radii, reflecting the differenc es in gas distributions in our model halos. Right: Comparison of the 
average Lyo: profile with observational results in|Momose et al.|(|2014|). Beyond ~ 40 kpc systematic effects in the image stacking analysis 
become important, indicated by the open circles, i lie shaded region gives an idea of the spread in the model profile, obtained by excluding 
the LAE with the faintest or the brightest extended profiles from the average. 


images of 3 16 2 : ^ 3.1 LAEs. According to [Momose] 


et al. (2014), the data at i? < 40 kpc are reliable, while 
at larger radii systematic effects in the stacking analysis 
become significant compared to the signal (see the top- 
middle panel in their Figure 8). We mark such a transi¬ 
tion by using filled circles at i? < 40 kpc and open circles 
at > 40 kpc for the data points. The shaded region 
around the mean model profile quantifies the uncertainty. 
The upper (lower) boundary is derived by excluding the 
LAE with the lowest (highest) surface brightness and av¬ 
eraging over the other eight LAEs. This serves to only 
provide some rough idea on the variation level of the 
mean profile, given the small number of model LAEs in 
our analysis. 

On small scales {R < 15 kpc), the model profile 
matches the observed profile extremely well, which is 
striking. At first glance one may attribute this to co¬ 
incidence, since we do not intend to fit the observed pro¬ 
file and we do not have any free parameters to adjust in 
our model. We directly use the Lya emissivity and gas 
distribution in the simulation. The only two changes we 
apply in the model besides the radiative transfer calcu¬ 
lation are “effective” dust and IGM absorption. The “ef¬ 
fective” Lya e xtinction optical dept h is the same as the 
one adopted in Cen & Zheng (2013) for studying LABs, 
which suppresses the initial intrinsic Lya emission. It 
aims at accounting for any uncertainties in the galaxy 
formation simulation. We tie it to the star formation 
rate and th e relation is fixed by considering halos above 
lO^^Mr:^ inICen & Zhengl (2013). We apply a mean ab- 

L\ 


sorption (scattering) tor Lya photons escaping the grid 
fro m the IGM outside of the simulation box. 

In Gen & Zheng| (|2013 ), the observed Lya luminosity- 
size relation of LABs and the Lya luminosity function 
of LABs are reproduced by our radiative transfer mod¬ 
eling. So it may not be too surprising that the similar 
model also provides a good match to the Lya emission in 
lower mass halos (10 ^^‘^Mq). Since the surface bright¬ 


ness profile at the central part is largely determined by 
the PSF, the agreement with the observation means that 
the central Lya luminosity in our model happens to be 
similar to the avera ge luminosity of the 31 6 LAEs in the 
stacking analysis in Momose et al. (2014), which by all 
means is an encouraging sign. 

At larger radii {R > 15 kpc), the model curve is slightly 
higher (at a factor of two level) than the observation. 
Given the small number of LAEs and the lack of adjust¬ 
ment in the model, the agreement to the observation still 
appears remarkable, in particular if the uncertainties in 
the data points and in the model curve are taken into 
account (keeping in mind that the data may suffer from 
significant systematic bias at i? > 40 kpc). Both the 
shape and extent of the LAH are reasonably reproduced. 

As a whole, our model mean Lya surface brightness 
profile, effectively computed from stacking Lya images 
of nine LAEs viewed along many different directions, 
shows good agreement with stacking analysis from ob¬ 
served ^ ^ 3.1 LAEs, from the central part to the diffuse 
LAH extended to ^ 60 kpc. However, the systematic 
effect in the data analysis makes the comparison beyond 
^ 40 kpc less reliable, and the apparent disagreement at 
R> 60 kpc is not significant. We proceed to investigate 
the contributions from various components to the surface 
brightness profile to gain more insights. 


3.2. Decomposing the Lya Surface Brightness Profile 

We record the initial position of each photon, which 
makes it possible to separate the contributions to the 
surface brightness profile from photons originating at dif¬ 
ferent places in our simulation. 

In each halo, there is a central LAE with strong Lya 
emission. We attribute Lya photons launched within 
10 kpc of the halo center to the central LAE. There are 
also a few regions in the halo with high Lya emissivity, 
which are small star-forming galaxies around the central 
LAE. In the Lya images shown in Figure such regions 
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R(kpc) R(kpc) 


Fig. 4.— Decomposition of the total Lya surface brightness profile into different components of Lya emission. Left: The decomposition 
of the total profile into contributions f rom the centra l LAE (red curve), other high emission knots (blue curve), and background regions 
(green curve). Observational data from|Momose et al.|(|2014|) is shown in purple. Systematic bias is important for > 40 kpc, indicated by 
the open circles. Note that at radii larger tJian ~ iU kpc, tJie profile from the central LAE is unable to account for observations, with the 
knot and background profiles playing dominant roles. Right: The decomposition of the total profile into contributions from star formation 
and cooling emission. The top panel shows profiles for each emission type, given in surface brightness units of erg s“^ cm“^ arcsec 
Star forming emission is shown in blue, and cooling emission is shown in red. The bottom panel shows the fractional contribution that 
each emission type makes to the total profile. 



r(kpc) 



Fig. 5. — Spatial distribution of the identified high emission knots in our simulation. Left: The distribution of 3D distances of the knots 
to the box center. The solid line shows a cumulative distribution function (CDF) of the 3D distances, with the dashed curve corresponding 
to a number density profile n{r) = 0 for r < 20 kpc and n{r) oc for r > 20 kpc. Right: The distribution of 2D projected distances of 
the knots to the box center. The solid line shows the CDF of the 2D distances and the dashed curve is for a number density distribution 
oc R~^. 


appear as relatively isolated peaks with surface bright¬ 
ness above ^ 10“^^erg s“^ cm“^ arcsec”^. The majority 


of them are below the det ection threshold for ty pical LAE 
surveys, as is the case in Momose et ah (2014). In three 
of the nine LAEs, a few of the high emissivity regions 
can reach the detection threshold, and would show up 
as isolated LAEs around the central LAEs. We refer to 
the high emissivity regions as “knots” and associate to 
each knot Lya photons launched within 10 kpc of its cen¬ 
ter. Photons that belong to neither the central LAE nor 


the knots are identified as being emitted from the back¬ 
ground of the simulation box. Most of them come from 
small clumps of gas that possess low rates of star forma¬ 
tion. Clearly the distinction between the knots and the 
background depends on our choice, which can be arbi¬ 
trary. However, our separation here serves the purpose of 
obtaining a rough idea on how Lya photons from differ¬ 
ent physical regions contribute to the surface brightness 
profile. 

The left panel of Eigure shows the decomposition 
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Fig. 6. — A further decomposition of the averaged surface brightness profile into contributions from star formation (blue curves) and 
cooling radiation (red curves). From left to right, the panels show profiles for photons produced within 10 kpc of the central LAEs, within 
the knot regions, and from the background areas. The top panels show surface brightness in units of erg s“^ cm“^ arcsec“^, while the 
bottom panels show the relative contributions. 


of the mean surface brightness profile (black) into con¬ 
tributions from Lyof photons originating in the central 
LAE (red), the knots regions (blue), and the background 
regions (green). 

After the radiative transfer, Lya photons originating 
from the central LAE appear to peak around the cen¬ 
tral region, following the PSE. The entire amplitude of 
the overall surface brightness profile at i? < 10 kpc 
comes from this component. The scatterings of photons 
with neutral hydrogen atoms in the circumgalactic and 
intergalactic media lead to an extended profile beyond 
^ 15 kpc. The profile drops toward large radii, roughly 
following which is too steep to account for the 

LAHs seen in both the model and observations. 

Lyof photons from the knots and background regions 
make comparable contributions (within a factor of about 
two) to the overall surface brightness profile at scales 
above ^ 20 kpc. They have similar profiles, which in turn 
are similar to that of the LAH and are flatter than the 
extended profile from scattered photons from the central 
LAE. Together, they dominate the profile at i? > 20 kpc. 

The above decomposition leads to the interesting im¬ 
plication that the observed extended Lya emission of 
LAHs is largely caused by emission from regions of low 
star formation rates spatially distributed inside the host 
dark matter halos of the central LAEs, and that pho¬ 
tons diffusing out from the central LAE as a result of 
the radiative transfer process play only a secondary role 
in producing the observed extended emission. 

We also decompose the surface brightness profile into 
contributions from Lya photons generated by star forma¬ 
tion and cooling radiation, as shown in the right panel 
of Eigure|^ The profiles from the two contributions are 
similar, but the photons from cooling radiation always 
sub-dominate, making about 30-40% of the total Lya 
light in the extended LAHs. 

Because the identified high-emission knots contribute 
significantly to the surface brightness profile at large 
radii, it is important to examine their spatial distribu¬ 


tion. The left panel of Eigure shows the cumulative 
distribution function (CDE) of the 3-dimensional (3D) 
distance r from the knots to the center for all the nine 
LAEs. Since we define the radius of each knot region 
and the central LAE to be 10 kpc, any high-emission ar¬ 
eas within 20 kpcof the center will not be identified as 
independent knots in our analysis. Therefore, the CDE 
curve starts at r = 20 kpc. The dashed curve corre¬ 
sponds to a number density distribution n(r) oc (for 
r > 20 kpc and 0 for r < 20 kpc). The plot shows that 
the high-emission knots closely follow a singular isother¬ 
mal distribution up to at least ^ 3i?vir- 

As we study the surface brightness profile, it would be 
more illustrative to examine the projected distribution 
of the knots. Eor this purpose, we chose a large number 
of isotropically distributed viewing directions. Eor each 
viewing direction and each LAE, we record the projected 
radius to the center for each knot. The right panel of Eig¬ 
ure shows the CDE for the 2D projected radii R. The 
dashed curve is the CDE for a surface number density 
that follows which appears to be a reasonable de¬ 

scription of the distribution of knots. This is consistent 
with the CDE of the 3D radii. The distribution of knots 
explains the shallow slope in the mean surface brightness 
profile seen in the model or observed LAHs (e.g., Eigure 
and Eigure®, which has a slope around —1. 

Eigure® snows the decomposition of our surface bright¬ 
ness profiles into contributions from star formation and 
cooling radiation for the central LAE, the identified knot 
regions, and the background. Eor Lya photons produced 
in the central LAE, cooling radiation makes up about 
1/3 of the observed Lya emission near the center, and 
its contribution drops to 15% at i? ^ 40 kpc. Eor Lya 
photons produced in the knot regions, cooling radiation 
and star formation contributions are comparable, and for 
those produced in background, the star formation con¬ 
tribution dominates. The latter two components depend 
on how we define knot regions. If we choose a lower 
threshold to define knots, star formation would become 











the dominant mechanism in producing Lya photons in 
knot regions. 

3.3. Possible Constraints from the UV Surfaee 
Brightness Profile 

The Lyof surface brightness profile from our model 
shows an encouraging agreement with the data. Besides 
Lya , the stacking analysis is als o performed for UV im¬ 
ages (e.g., Momose et al. 2014). UV photons are pro¬ 
duced froni star formatiom Unlike Lya photons, they 
do not interact with neutral hydrogen through resonant 
scattering. Instead, they escape directly from their point 
of creation (modulated by dust extinction), which allows 
them to serve as a tool to map out regions of star for¬ 
mation. The UV profile of LAEs can therefore provide 
complementary information about LAHs and can be used 
to further constrain the origin of LAHs. 

We convert the SFR to UV luminosity (at rest- 
frame 1500A) using the prescri ption Luv = 8 x 

10 ^^[SFR/(M 0 yr“^)]ergs“^Hz“^ (Madau et al.||l998|). 


It assumes a Salpeter initial mass function (IIVH') of stars 
and solar metallicity. 

While dust extinction should also be included to obtain 
the observed UV luminosity, it is degenerated with the 
above assumption about the stellar population and to a 
less degree with the star formation history. For example, 
a sub-solar metallicity (or a different IMF, e.g., Chabrier 
IMF) results in a higher SFR-to-UV conversion factor, 
and the combined effect of the metallicity and IMF can 
lead to a factor of a few increase in the conversio n (e.g., 
Leitherer et al. 1999 Madau & Dickinson 2014). The 
dust extinction works in the direction to bring tne con¬ 
version factor toward the above value. Given such un¬ 
certainties in the model, we can simply adopt the above 
conversion to proceed, and rescale the model UV profile 
if necessary to fit the observed profile. 

We create UV photons based on the SFR in each cell 
and assign a random escape direction for each photon. 
Figure [^displays the UV images of the nine model galax¬ 
ies observed along the same direction as in Figure To 
study the average UV light distribution seen in the st ack- 
ing analysis, we follow the method described in § | 2 . 2 | 
to produce the average UV surface brightness profile for 
each model LAE. 

The black curve in the left panel of Figure shows 
the UV profile predicted by our model. At small radii 
{R < 15 kpc), the model curve is almost right on top of 
the observed profile. 

Such a coincidence indicates that our SFR-to-UV con¬ 
version factor is about right in reflecting a combined 
effect of stellar population, metallicity, extinction, and 
model SFR, even though there are uncertainties in each 
component and the overall model is approximate. As 
such, we make no adjustments to our initial conversion 
factor. 

The model reproduces the central UV luminosity, and 
the shape of the profile simply follows that of the PSF. 

At large radii {R > 15 kpc), our model shows an ex¬ 
tended UV halo (left panel of Figure 1^. This is not 
unexpected, given that emission product from star for¬ 
mation in the outer halo makes a substantial contribution 
to the extended Lya profile (Figure]^. The model UV 
profile is in apparent tension with observations (Momose 


|et al.||2014|), where little evidence is shown for such ex- 
tended UV halos (see the data points in the left panel of 
Figure [^. 

We note that there are some residuals of sky subtrac¬ 
tion fo und in the composite UV images of Momose et al. 
(2014). To quantify the significance of the apparent ten- 
sion between the model and observation, we evaluate 
the sky subtraction systematics in the average UV sur¬ 
face brightness profile. We find that there is a signal 
of sky over-subtraction at the surface brightness level of 
3.0 X 10“^^ergs“^ cm“^ arcsec“^, corresponding to 
a maximum negative value of the UV surface brightness 
estimate of the LAH (Momose et al. 2015, in prepara¬ 
tion). Therefore, we conclude that the average UV profile 
below the above level is subject to the influence of the 
systematics. Note that the effect of sky over-subtraction 
systematics is canceled out in the Lya surface brightness 
profile of the LAH, because the Lya profile is obtained 
by taking the difference of the composite broadband and 
narrowband images, which have the same level of sky 
over-subtraction. 

The dashed line in the left panel of Figure shows 
the sky-subtraction systematic effect at the surface 
brightness level of 3.0x 10“^^erg s“^ cm“^ Hz“^ arcsec“^, 
which can be regarded as an upper bound for extended 
UV profile. Our model UV profile (black curve) appears 
to be at a similar level. Therefore, by accounting for the 
systematics in the data, no significant tension is found 
between the model and the observation. Given the sky- 


subtraction system atics, the UV profile obtained in Mo 


mose e t al.| ( 2Q14[ ) (data points in the left panel of our 


h'igure|^ represents a lower bound, which essentially fol¬ 
lows the PSF and is determined by the central star for¬ 
mation. To match this case, we construct a modified 
model by removing star formation in the outer halo re¬ 
gions in our model. It is not surprising to see that the 
modified model profile produced in this way (blue curve 
in the left panel of Figure]^ follows the PSF, as with the 
data points. 

The default model profile (black curve) and the modi¬ 
fied one (blue curve) should be able to bracket all possi¬ 
ble cases of the extended UV profile. Improved measure¬ 
ments of the UV profile with well-controlled systematics 
would provide important information on the amount of 
star formation in the outer halo. This in turn would 
improve our understanding of the origin of the extended 
Lya halo, along with its partition into contributions from 
star formation and cooling radiation. With the current 
situation, the above two boundary cases allow us to infer 
the range of the relative contributions of star formation 
and cooling radiation to the Lya surface brightness pro¬ 
file. 

The black curve in the right panel of Figure corre¬ 
sponds to the Lya profile from our default model. We 
find that the total Lya luminosity within the projected 
radius R < Rvir = 56 kpc can be broken down into the 
following contributions: 33% from star forming photons 
produced in the central galaxy, 15% from cooling radia¬ 
tion emitted from the central galaxy (r < 10 kpc), 28% 
from star forming photons produced in the outer halo, 
and 25% from cooling radiation in the outer halo. In the 
extended part of the profile (defined as emission observed 
at 15 kpc ^ R ^ 56 kpc), the fraction of photons from 
cooling radiation is about 42%. 
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Fig. 7.— UV surface brightness images for all nine model LAEs in our original analysis. Each image is 224 kpc (physical) on a side and 
has been smoothed by a 2D Gaussian kernel with a FWHM of 1.32" to match observations. The dashed circle in e ach panel has a radiu s 
of 40 kpc, roughly corresponding to the radius that systematic effects become important in the stacking analysis in|Momose et al.| (|2014|). 
The observed profile can potentially be used to put constraints on the clustered UV sources around the central galaxies, as discussed in 
the text. 


The Lya profile after removing star formation from 
the outer halo is shown as the blue curve in the right 
panel of Figure The effect is not drastic, and the Lya 
profile for this case drops to a level in even better agree¬ 
ment with observations. For this modified profile, the 
average component contributions to the total Lya lumi¬ 
nosity within projected radius R < i?vir = 56 kpc are 
45% from star forming photons produced in the central 
galaxy, 20% from cooling radiation within and around 
the central galaxy (r < 10 kpc), and 35% from cooling 
radiation in the outer halo. On average, cooling radia¬ 
tion can now contribute to about half of the total Lya 
luminosity. If we focus on the extended part of the profile 
(15 kpc ^ ^ 56 kpc), the fraction of the cooing radia¬ 

tion contribution is about 75%. With this prescription, 
the extended Lya halo is dominated by cooling radiation. 

Although the best currently available data is limited by 
systematics, the UV profile in combination with the Lya 
profile can help to constrain the nature of LAHs. Our 
investigation implies that cooling radiation in the outer 
halo may play a significant role in forming extended Lya 
halos (e.g., contributing more than half of the emission), 
a prediction that can be tested with tighter observational 
constraints on the UV profile. The caveats and more 
discussions are presented in the next section. 

4. SUMMARY AND DISCUSSION 

We perform Lya radiative transfer modeling of 2 ; = 3.1 
LAEs with a high resolution hydrodynamic cosmological 


galaxy formation simulation to study the extended Lya- 
emitting halos recently discovered in observation from 
image stacking analysis. We develop a method to com¬ 
pute the mean surface brightness profile from averaging 
over many different viewing directions. We consider nine 
model LAEs residing in halos of and find their 

mean Lya surface brightness profile to b e in remarkable 
agreem ent with the observed profile in |Momose et al.| 
(2014), at both the central and extended parts. 

T'o investigate the origin of the extended Lya emis¬ 
sion, we decompose the profile into contributions from 
Lya photons produced in different regions, which include 
the central LAE in each simulation box, dense regions of 
high star formation activity spatially separated from the 
central LAE (dubbed as “knots”), and faint background 
areas. The latter two outer halo components are associ¬ 
ated with satellite galaxies or tidally stripped materials 
in the halo. 

Lya photons originating near the halo center (from 
both star formation and cooling radiation) but scattered 
to large radii by the hydrogen atoms in the CG M do pro¬ 


duce an extended Lya halo , as pr edicted by, e.g. , Laursen 
& Sommer-Larsen ( 2QQ7[ lp and Zheng et al. ( |2Qllb ) 


However, our radiative transfer model with the high- 


^ Cooling radiation in t he outer halo is also included in |Laursen| 
|& Sommer-Larsen| (|2007|). However, there is no discussion on its 
contribution to the Lam for an obvious reason - LAHs had not yet 
been discovered in observations. 
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Fig. 8.— Average UV and Lyg surface b rightness profiles and the effect of star formation in the outer halo. In each panel, observational 
data are taken from |Momose et al.| (|2014|) with open circles representing the region where the data starts to be limited by systematics. 
Left: Comparison of our U V prohles with observational data. The black curve is the average profile from our initial model. The blue curve 
shows the effect of removing star formation in the outer halo, keeping only the star forming emission from the central LAE. The dashed line 
is drawn at 3.0 x 10“^^ergs“^ cm“^ Hz“^ arcsec“^, which is the level of systematic effects in the image stacking analysis (see the text). 
This line represents an upper limit for any extended UV emission that actually occurs. Right: The corresponding changes in the model 
Lya profile. The black curve shows our total Lya profile, and the blue curve shows the effect of removing star formation emission from the 
outer halo. The gray shaded region in the right panel gives an idea of the spread of the profile, obtained by excluding the LAE with the 
faintest or the brightest extended profiles from the average. 

resolution galaxy formation simulation shows that such 
a contribution alone is not able to explain the surface 
brightness level of the observed LAH (e.g., lower by a 
factor of 10 around R = 40 kpc). Instead, the extended 
LAH is dominated by emission from the knots and the 
background regions of the outer halo. The result im¬ 
plies that scattering of Lya photons from bright, central 
sources is less important in forming LAHs than previ¬ 
ously thought. 

Certainly the exact profile created by the scattered, 
centrally produced photons should depend on the density 
and velocity distribution of the circumgalactic gas. Some 
analytic models with clumpy CGM and decelerating out¬ 
flows can produce scattered Lya halos at the observed 
surface brightne ss level (e.g., Steidel et al.|2011 Dijkstra 
fc Kramer| 2Q12 ). It is worth investigating the contribu¬ 


tion from such scattered halos with more realistic CGM 
distributions from high-resolution galaxy formation sim- 


feedback (e.^ 

., Suresh et al. 2015 Muratov et al.||20I5i 

Kimm et al. 

2 OI 5 I). By performing a test with a lower- 


we also find that the scattered profile shows a weak de¬ 
pendence on the resolution, becoming slightly more ex¬ 
tended with higher resolution. However, the scattered 
profile still drops rapidly toward large radii, remaining 
unlikely to account for the observed one. 

Lya emission from the outer halo has two components 
- gravitational cooling radiation and emission from star 
formation. In the simulation, we find that star forma¬ 
tion slightly dominates, and cooling radiation makes a 
substantial contribution. Our model predicts the exis¬ 
tence of an extended UV halo at a brightness level of 
^ 3 X I0“^^ergs“^ cm“^ Hz“^ arcsec“^, which is right at 
the limit of the sky-subtraction systematics in observa¬ 


tional data (Momose et al. 2014). If actual UV halos 
(from improved data analysis) are significantly dimmer 
than this, we will need to investigate how to suppress the 
UV profile in the model (e.g., identifying the likely cause 
in the simulation for the over-prediction of star forma¬ 
tion in the outer halo or studying the dust attenuation 
effect in the outer halo). Although the best currently 
available UV profile measurement does not serve as a ro¬ 
bust constraint as a result of systematics, we find that 
the extended Lya profile becomes even better reproduced 
after removing star formation in the outer halo in order 
to suppress the UV profile. 

Taken at face value, our investigation shows that our 
initial model can ex plain both the observed Lya and UV 
profiles of LAEs in Momose et al. (2014), from small 
radii (< 15 kpc) to large radii (up to ^ 80 kpc). The 
agreements between the model and data are excellent. 
This is remarkable, especially given that we do not in¬ 
tend to fit the profiles by tuning parameters. Depending 
on the accuracy of our star-forming recipe and the loose 
constraints from the UV halo, we find that cooling radi¬ 
ation can contribute 40-55% of the total Lya luminosity 
within projected radius R < i?vir, where i?vir = 56 kpc 
(^ 7.2") is the virial radius of the LAE host halo. Eor the 
diffuse LAH, which is usually buried in sky noise for in¬ 
dividual LAEs, the contribution from cooling radiation 
is more substantial, making up about 42-75% (within 
15 kpc ^ ^ 56 kpc). 

Gravitational cooling radiation from accretion of gas 
is a process expected to occur during galaxy formation, 
mainly in the form of Lya emission from collisional exci¬ 
tation and de-e xcitation of hydrog en atoms in gas around 


2 X 10^ K (e.g., Eardal et al.|2QQI ). Many previous stud¬ 


ies of cooling radiation with analytic calculations and hy- 
drodynamical simulations focus on investigating it as a 
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possible mechanism to explain Lyg blobs (e.g., Haiman 
et al.||2QQQ| [Fardal et al. |2QQ1| [Furlanetto et al.||2UU^ 

L( 


Yang et 'nl|2QQ6[ IDijkstra fc Loeb||2QQ9[ [Goerdt et ^ 

J Faucher-Giguere et al.||2010|), which are more lu- 




minoiis~than LAH s (b teidel et al.|2QQQ). As shown by 
Yang et al. (2006) ana b'aucher-Giguere et ah (2010), an 
accurate prediction of the cooling Lya emission relies on 
an accurate treatment of the self-shielding effect for the 
ionizing photons, which affects the ionization and ther¬ 
mal states of the accreted gas. In the simulation used 
in this work, self-shielding correction is performed on- 
the-fly. There could be small variations in the predicted 
cooling radiation with different self-shielding correction 
methods, but it is unlikely to remove the cooling radia¬ 
tion signal in our model, which is significant regardless 
of the accuracy of our star-forming recipe. 

One caveat to keep in mind is that our results in this 
paper are based on radiative transfer modeling of nine 
simulated galaxies in halos of IO^^’^Mq. First, our anal¬ 
ysis suffers from small number statistics. While we at¬ 
tempt to make full use of the nine galaxies by obtaining 
the mean surface brightness profile from averaging all 
viewing angles (effectively creating a much larger sam¬ 
ple for stacking), modeling more galaxies definitely helps 
the study of LAHs. More galaxies are also needed to ex¬ 
plore the dependence of LAHs on the properties of galax - 
ies and their environments (e.g., Matsuda et al. 2012). 
Second, the mass of halos (IO^^’^Mq) considered in this 
work seems to be on the high end of LAE-hosting halos. 
The LAE halo masses in ferred from cluster ing analysis 
are typically (Ouchi et al. 2010). Clearly it 

is necessary to investigate how the Lya and UV surface 
brightness profiles and their decomposition into cooling 


and star forming components vary with halo mass, and 
to make comparisons with data especially as better ob¬ 
servational constraints on the UV profile become avail- 
able. As an example of th e potential impact of halo mass. 


Rosdahl & Blaizot (2012) find that the extent of cooling 
radiation in the outer halo is dependent on halo mass, re¬ 
sulting from a positive correlation between the efficiency 
of cold streaming accretion and halo mass. Additional 
radiative transfer modeling of galaxies in lower mass ha¬ 
los comparing emission from the central LAE with star 
formation and cooling radiation in the outer halo will 
elucidate to what extent our results in this paper hold. 

We plan to carry out studies related to the halo mass 
and redshift evolution of LAHs in order to obtain a better 
understanding of their origins and to learn more about 
the CGM and galaxy formation. 


We thank the anonymous referee for constructive 
comments and Mark Dijkstra for helpful discussions. 
E.L., Z.Z., and R.S. are supported by NASA grant 
NNX14AC89G and NSE grant AST-1208891. The sup¬ 
port and resources from the Genter for High Perfor¬ 
mance Gomputing at the University of Utah are grate¬ 
fully acknowledged. Gomputing resources were in part 
provided by the NASA High-End Gomputing (HEC) 
Program through the NASA Advanced Supercomput- 
ing (NAS) Division at Ames Research Genter. R.G. 
is supported in part by grants NASA NNX11AI23G 
and NNX12AE91G. This work was supported by World 
Premier International Research Genter Initiative (WPI 
Initiative), MEXT, Japan, and KAKENHI (23244025) 
Grantin-Aid for Scientific Research (A) through Japan 
Society for the Promotion of Science (JSPS). 


REFERENCES 


Barnes, L. A., & Haehnelt, M. G. 2010, MNRAS, 403, 870 
Barnes, L. A., Haehnelt, M. G., Tescari, E., & Viel, M. 2011, 
MNRAS, 416, 1723 

Brinchmann, J., Chariot, S., White, S. D. M., et al. 2004, 
MNRAS, 351, 1151 

Bryan, G. L., & Norman, M. L. 2000, Institute for Mathematics 
and Its Applications, 117, 165 
Gen, R. 2011, ApJ, 742, L33 
—. 2012, ApJ, 748, 121 
Gen, R., & Zheng, Z. 2013, ApJ, 775, 112 
Dijkstra, M., & Kramer, R. 2012, MNRAS, 424, 1672 
Dijkstra, M., & Loeb, A. 2009, MNRAS, 400, 1109 
Fardal, M. A., Katz, N., Gardner, J. P., et al. 2001, ApJ, 562, 605 
Faucher-Giguere, C.-A., Keres, D., Dijkstra, M., Hernquist, L., & 
Zaldarriaga, M. 2010, ApJ, 725, 633 
Feldmeier, J. J., Hagen, A., Ciardullo, R., et al. 2013, ApJ, 776, 
75 

Furlanetto, S. R., Schaye, J., Springel, V., & Hernquist, L. 2005, 
ApJ, 622, 7 

Fynbo, J. P. U., Ledoux, C., Moller, P., Thomsen, B., & Burud, I. 
2003, A&A, 407, 147 

Gawiser, E., Francke, H., Lai, K., et al. 2007, ApJ, 671, 278 
Goerdt, T., Dekel, A., Sternberg, A., et al. 2010, MNRAS, 407, 
613 

Guaita, L., Gawiser, E., Padilla, N., et al. 2010, ApJ, 714, 255 
Haiman, Z., Spaans, M., & Quataert, E. 2000, ApJ, 537, L5 
Jeeson-Daniel, A., Ciardi, B., Maio, U., et al. 2012, MNRAS, 424, 
2193 

Jiang, L., Egami, E., Fan, X., et al. 2013, ApJ, 773, 153 
Joung, M. R., Cen, R., & Bryan, G. L. 2009, ApJ, 692, LI 
Kimm, T., Cen, R., Devriendt, J., Dubois, Y., & Slyz, A. 2015, 
ArXiv e-prints, arXiv: 1501.05655 
Laursen, P., & Sommer-Larsen, J. 2007, ApJ, 657, L69 


Laursen, P., Sommer-Larsen, J., Sz Andersen, A. C. 2009, ApJ, 
704, 1640 

Leitherer, C., Schaerer, D., Goldader, J. D., et al. 1999, ApJS, 
123, 3 

Madau, P., Pozzetti, L., & Dickinson, M. 1998, ApJ, 498, 106 
Madau, P., & Dickinson, M. 2014, ARA&A, 52, 415 
Matsuda, Y., Yamada, T., Hayashino, T., et al. 2012, MNRAS, 
425, 878 

Momose, R., Ouchi, M., Nakajima, K., et al. 2014, MNRAS, 442, 
110 

Muratov, A. L., Keres, D., Faucher-Giguere, C.-A., et al. 2015, 
ArXiv e-prints, arXiv:1501.03155 
Ouchi, M., Shimasaku, K., Akiyama, M., et al. 2008, ApJS, 176, 
301 

Ouchi, M., Shimasaku, K., Furusawa, H., et al. 2010, ApJ, 723, 
869 

Partridge, R. B., & Peebles, P. J. E. 1967, ApJ, 147, 868 
Rhoads, J. E., Dey, A., Malhotra, S., et al. 2003, AJ, 125, 1006 
Rosdahl, J., & Blaizot, J. 2012, MNRAS, 423, 344 
Steidel, C. C., Adelberger, K. L., Shapley, A. E., et al. 2000, ApJ, 
532, 170 

Steidel, C. C., Bogosavljevic, M., Shapley, A. E., et al. 2011, ApJ, 
736, 160 

Suresh, J., Bird, S., Vogelsberger, M., et al. 2015, MNRAS, 448, 
895 

Tasitsiomi, A. 2006, ApJ, 645, 792 

Yang, Y., Zabludoff, A. L, Dave, R., et al. 2006, ApJ, 640, 539 
Zahid, H. J., Yates, R. M., Kewley, L. J., Sz Kudritzki, R. P. 

2013, ApJ, 763, 92 

Zheng, Z., Cen, R., Trac, H., & Miralda-Escude, J. 2010, ApJ, 
716, 574 

—. 2011a, ApJ, 726, 38 

Zheng, Z., Cen, R., Weinberg, D., Trac, H., & Miralda-Escude, J. 
2011b, ApJ, 739, 62 


















































12 


Zheng, Z., & Miralda-Escude, J. 2002, ApJ, 578, 33 
Zheng, Z., & Wallace, J. 2014, ApJ, 794, 116 



